

********************************************************************************
************************************ Figure 1 **********************************
********************************************************************************
use data_1, clear

reg support_milforce i.war_treatment_4cat, r     
margins i.war_treatment_4cat, level(95) asobserved 
marginsplot, ///
plotopts(mcolor(black) lpattern(blank) lcolor(black) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ciopts(color(black) lcolor(black) lwidth(medthick)) recastci(rspike) ///
level(95) graphregion(fcolor(white) lcolor(white) margin(medium)) ///
plotr(m(medsmall)) ///
ylabel(0(1)4, nogrid labsize(medsmall))  ///
xlabel(-0.5 " " 0 "None" 1 "Prov" 2 "Prov+De" 3 "Prov+Esc" 3.5 " ", nogrid labsize(*1.2) noticks) ///
xtitle(" " , margin(0 0 0 0) color(white)) title(" ", color(black)) ytitle(" ",  margin(0 0 0 0) color(black)) xsize(15) ysize(20) ///
legend(region(lstyle(none) lcolor(white)) order(2 "All") size(*1.2) symxsize(*.4) symysize(*.4) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
addplot( ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==0,  discrete width(1) freq yaxis(3) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)900, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)900, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==2,  discrete width(1) freq yaxis(3) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)900, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==3,  discrete width(1) freq yaxis(3) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)900, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==4,  discrete width(1) freq yaxis(3) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)900, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==0,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==2,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==3,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==4,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==0,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==2,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==3,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==4,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==0,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==2,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==3,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==4,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(*1.2) nogrid) ylabel(0(100)1800, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(3)) ytitle("", size(zero) margin(zero) axis(3))) ///
, ///
legend(region(lstyle(none) lcolor(white)) order(2 "All") size(*1) symxsize(*.3) symysize(*.3) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
)

graph save all, replace


reg support_milforce i.war_treatment_4cat##c.putin_approve_dum_c gender i.profile_age1 i.Qregion_federal_districts_recode i.income_recoded, r     
margins, at(war_treatment_4cat = (0 1 2 3) putin_approve_dum_c = (0 1)) level(95) asobserved 
marginsplot, ///
plotopts(mcolor(eltblue) lpattern(blank) lcolor(eltblue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ciopts(color(eltblue) lcolor(eltblue) lwidth(medthick)) recastci(rspike) ///
plot2opts(mcolor(blue) lpattern(blank) lcolor(blue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ci2opts(color(blue) lcolor(blue) lwidth(medthick))  ///
level(95) graphregion(fcolor(white) lcolor(white) margin(medium)) ///
ylabel(0(1)4, nogrid labsize(medsmall))  ///
xlabel(-0.5 " " 0 "None" 1 "Prov" 2 "Prov+De" 3 "Prov+Esc" 3.5 " ", nogrid labsize(*1.2) noticks) ///
xtitle(" " , margin(0 0 0 0) color(white)) title(" ", color(black)) ytitle(" ", margin(0 0 0 0) color(white)) xsize(15) ysize(20) ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1.2) symxsize(*.4) symysize(*.4) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
addplot( ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ylabel(0(1)4, nogrid labsize(*1.2))  xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(3)) ytitle("", size(zero) margin(zero) axis(3))) ///
, ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1) symxsize(*.3) symysize(*.3) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
)

graph save sup, replace


graph combine all.gph sup.gph, xsize(20) ysize(10) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) b1("Treatments" , size(*1.2) margin(0 0 0 0) color(black)) l1("Support for war" , size(*1.2) margin(0 0 0 0) color(black)) rows(1) cols(3) 

graph export figure_1.png, width(4000) replace
graph export figure_1_pdf.pdf, replace



********************************************************************************
************************************ Figure 2 **********************************
********************************************************************************
use data_1, clear

reg support_milforce i.war_treatment_4cat##c.putin_approve_dum_c gender i.profile_age1 i.Qregion_federal_districts_recode i.income_recoded if treatment_culture==1 | war_treatment==0, r     
margins, at(war_treatment_4cat = (0 1 2 3) putin_approve_dum_c = (0 1)) level(95) asobserved 
marginsplot, ///
plotopts(mcolor(eltblue) lpattern(blank) lcolor(eltblue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ciopts(color(eltblue) lcolor(eltblue) lwidth(medthick)) recastci(rspike) ///
plot2opts(mcolor(blue) lpattern(blank) lcolor(blue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ci2opts(color(blue) lcolor(blue) lwidth(medthick))  ///
level(95) graphregion(fcolor(white) lcolor(white) margin(medium)) ///
ylabel(0(1)4, nogrid labsize(medsmall))  ///
xlabel(-0.5 " " 0 "None" 1 "Prov" 2 "Prov+De" 3 "Prov+Esc" 3.5 " ", nogrid labsize(*1.2) noticks) ///
xtitle(" " , margin(0 0 0 0) color(white)) title(" ", color(black)) ytitle(" ", margin(0 0 0 0) color(white)) xsize(15) ysize(20) ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1.2) symxsize(*.4) symysize(*.4) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
addplot( ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==0 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==1 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==2 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==3 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==4 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==0 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==1 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==2 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==3 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==4 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==0 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==1 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==2 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==3 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==4 & putin_approve_dum_c==1 & treatment_culture==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, nogrid labsize(*1.2))  xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("Culture treatments", size(vlarge) margin(0 0 0 0) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(3)) ytitle("", size(zero) margin(zero) axis(3))) ///
, ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1) symxsize(*.3) symysize(*.3) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
)

graph save cul_put, replace


reg support_milforce i.war_treatment_4cat##c.putin_approve_dum_c gender i.profile_age1 i.Qregion_federal_districts_recode i.income_recoded if treatment_military==1 | war_treatment==0, r     
margins, at(war_treatment_4cat = (0 1 2 3) putin_approve_dum_c = (0 1)) level(95) asobserved 
marginsplot, ///
plotopts(mcolor(eltblue) lpattern(blank) lcolor(eltblue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ciopts(color(eltblue) lcolor(eltblue) lwidth(medthick)) recastci(rspike) ///
plot2opts(mcolor(blue) lpattern(blank) lcolor(blue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ci2opts(color(blue) lcolor(blue) lwidth(medthick))  ///
level(95) graphregion(fcolor(white) lcolor(white) margin(medium)) ///
ylabel(0(1)4, nogrid labsize(medsmall))  ///
xlabel(-0.5 " " 0 "None" 1 "Prov" 2 "Prov+De" 3 "Prov+Esc" 3.5 " ", nogrid labsize(*1.2) noticks) ///
xtitle(" " , margin(0 0 0 0) color(white)) title(" ", color(black)) ytitle(" ", margin(0 0 0 0) color(white)) xsize(15) ysize(20) ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1.2) symxsize(*.4) symysize(*.4) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
addplot( ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==0 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==1 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==2 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==3 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==4 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==0 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==1 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==2 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==3 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==4 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==0 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==1 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==2 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==3 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==4 & putin_approve_dum_c==1 & treatment_military==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, nogrid labsize(*1.2))  xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("Security treatments", size(vlarge) margin(0 0 0 0) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(3)) ytitle("", size(zero) margin(zero) axis(3))) ///
, ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1) symxsize(*.3) symysize(*.3) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
)


graph save mil_put, replace

graph combine  cul_put.gph mil_put.gph, xsize(20) ysize(10) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) l1("Support for war" , size(*1.2) margin(0 0 0 0) color(black)) rows(1) cols(3) 

graph export figure_2.png, width(4000) replace
graph export figure_2_pdf.pdf, replace


********************************************************************************
************************************ Figure 3 **********************************
********************************************************************************
use data_1, clear

reg support_milforce i.war_treatment_4cat##c.putin_approve_dum_c gender i.profile_age1 i.Qregion_federal_districts_recode i.income_recoded if treatment_military==1 | war_treatment==0, r     

graph pie if e(sample) & war_treatment_4cat==0, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))   ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("{bf:All}", size(vlarge) margin(0 4 0 0))  plotregion(margin(7.5 0 0 0)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white) margin(zero)) 
graph save none_all, replace

graph pie if e(sample) & war_treatment_4cat==1, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))    ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise plotregion(margin(9 0 0 0)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white) margin(zero)) 
graph save prov_all, replace

graph pie if e(sample) & war_treatment_4cat==2, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))   ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise  plotregion(margin(9 0 0 0)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white) margin(zero)) 
graph save de_all, replace

graph pie if e(sample) & war_treatment_4cat==3, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8)) ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise plotregion(margin(9 0 0 0)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white) margin(zero)) 
graph save es_all, replace


graph pie if e(sample) & putin_approve_dum_c==1 & war_treatment_4cat==0, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))   ///
plabel(_all percent, format(%4.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("{bf:Putin supporters}", size(vlarge) margin(0 4 0 0)) b1title("xxx", size(vlarge) margin(0 0 0 2) color(white)) plotregion(margin(medium)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save none_put, replace

graph pie if e(sample) & putin_approve_dum_c==1 & war_treatment_4cat==1, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))    ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("xxx", size(vlarge) margin(0 4 0 0) color(white)) b1title("xxx", size(vlarge) margin(0 0 0 2) color(white)) plotregion(margin(medium)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save prov_put, replace

graph pie if e(sample) & putin_approve_dum_c==1 & war_treatment_4cat==2, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))   ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("xxx", size(vlarge) margin(0 4 0 0) color(white)) b1title("xxx", size(vlarge) margin(0 0 0 2) color(white)) plotregion(margin(medium)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save de_put, replace

graph pie if e(sample) & putin_approve_dum_c==1 & war_treatment_4cat==3, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8)) ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("xxx", size(vlarge) margin(0 4 0 0) color(white)) b1title("xxx", size(vlarge) margin(0 0 0 2) color(white)) plotregion(margin(medium)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save es_put, replace


graph pie if e(sample) & putin_approve_dum_c==0 & war_treatment_4cat==0, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))    ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("{bf:Putin opponents}", size(vlarge) margin(0 4 0 0)) b1title("{bf:None}", size(vlarge) margin(0 0 0 2)) plotregion(margin(medium))  ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save none_opp, replace

graph pie if e(sample) & putin_approve_dum_c==0 & war_treatment_4cat==1, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))    ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("xxx", size(vlarge) margin(0 4 0 0) color(white)) b1title("{bf:Prov}", size(vlarge) margin(0 0 0 2))  plotregion(margin(medium))  ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save prov_opp, replace

graph pie if e(sample) & putin_approve_dum_c==0 & war_treatment_4cat==2, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))   ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("xxx", size(vlarge) margin(0 4 0 0) color(white)) b1title("{bf:Prov+De}", size(vlarge) margin(0 0 0 2))  plotregion(margin(medium)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save de_opp, replace

graph pie if e(sample) & putin_approve_dum_c==0 & war_treatment_4cat==3, over(support_milforce_3cat) pie(1, color(forest_green%60)) pie(2, color(orange%60)) pie(3, color(red%60) explode(8))   ///
plabel(_all percent, format(%5.0f) gap(-8.0) size(large) color(black)) title("",) noclockwise l1title("xxx", size(vlarge) margin(0 4 0 0) color(white)) b1title("{bf:Prov+Esc}", size(vlarge) margin(0 0 0 2))  plotregion(margin(medium)) ///
line(lcolor(black) lwidth(vthin)) intensity(inten70) legend(off) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) 
graph save es_opp, replace


graph combine none_all.gph prov_all.gph de_all.gph es_all.gph none_put.gph prov_put.gph de_put.gph es_put.gph none_opp.gph prov_opp.gph de_opp.gph es_opp.gph, xsize(20) ysize(10) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white) margin(zero)) plotregion(margin(zero)) rows(3) cols(4) imargin(zero)

graph export figure_3.png, width(4000) replace
graph export figure_3_pdf.pdf, replace



********************************************************************************
************************************ Figure 4 **********************************
********************************************************************************
use data_1, clear

reg support_milforce i.war_treatment_4cat##c.putin_approve_dum_c gender i.profile_age1 i.Qregion_federal_districts_recode i.income_recoded if treatment_country==1, r     
margins, at(war_treatment_4cat = (0 1 2 3) putin_approve_dum_c = (0 1)) level(95) asobserved 
marginsplot, ///
plotopts(mcolor(eltblue) lpattern(blank) lcolor(eltblue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ciopts(color(eltblue) lcolor(eltblue) lwidth(medthick)) recastci(rspike) ///
plot2opts(mcolor(blue) lpattern(blank) lcolor(blue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ci2opts(color(blue) lcolor(blue) lwidth(medthick))  ///
level(95) graphregion(fcolor(white) lcolor(white) margin(medium)) ///
ylabel(0(1)4, nogrid labsize(medsmall))  ///
xlabel(-0.5 " " 0 "None" 1 "Prov" 2 "Prov+De" 3 "Prov+Esc" 3.5 " ", nogrid labsize(*1.2) noticks) ///
xtitle(" " , margin(0 0 0 0) color(white)) title(" ", color(black)) ytitle(" ", margin(0 0 0 0) color(white)) xsize(15) ysize(20) ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1.2) symxsize(*.4) symysize(*.4) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
addplot( ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ylabel(0(1)4, nogrid labsize(*1.2))  xtitle("Latvia", size(vlarge) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(3)) ytitle("", size(zero) margin(zero) axis(3))) ///
, ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1) symxsize(*.3) symysize(*.3) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
)

graph save fig5_lat, replace


reg support_milforce i.war_treatment_4cat##c.putin_approve_dum_c gender i.profile_age1 i.Qregion_federal_districts_recode i.income_recoded if treatment_country==2, r     
margins, at(war_treatment_4cat = (0 1 2 3) putin_approve_dum_c = (0 1)) level(95) asobserved 
marginsplot, ///
plotopts(mcolor(eltblue) lpattern(blank) lcolor(eltblue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ciopts(color(eltblue) lcolor(eltblue) lwidth(medthick)) recastci(rspike) ///
plot2opts(mcolor(blue) lpattern(blank) lcolor(blue) lwidth(thin) msymbol(circle) msize(*0.75) mlwidth(*0.75)) ci2opts(color(blue) lcolor(blue) lwidth(medthick))  ///
level(95) graphregion(fcolor(white) lcolor(white) margin(medium)) ///
ylabel(0(1)4, nogrid labsize(medsmall))  ///
xlabel(-0.5 " " 0 "None" 1 "Prov" 2 "Prov+De" 3 "Prov+Esc" 3.5 " ", nogrid labsize(*1.2) noticks) ///
xtitle(" " , margin(0 0 0 0) color(white)) title(" ", color(black)) ytitle(" ", margin(0 0 0 0) color(white)) xsize(15) ysize(20) ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1.2) symxsize(*.4) symysize(*.4) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
addplot( ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t0 if e(sample) & war_treatment_4cat==0 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(3) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)346, labsize(vsmall) nolabels noticks axis(3))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t1 if e(sample) & war_treatment_4cat==1 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t2 if e(sample) & war_treatment_4cat==2 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==0 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sand%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==1 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(sandb%75) fintensity(75) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==2 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==3 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(orange_red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) xtitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2))) ///
(histogram support_milforce_t3 if e(sample) & war_treatment_4cat==3 & support_milforce==4 & putin_approve_dum_c==1,  discrete width(1) freq yaxis(2) xaxis(2) color(red%60) fintensity(60) xlabel(-5(1)35, labsize(vsmall) nolabels noticks axis(2)) ylabel(0(1)4, labsize(medsmall) nogrid) ylabel(0(10)692, labsize(vsmall) nolabels noticks axis(2)) xscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ylabel(0(1)4, nogrid labsize(*1.2))  xtitle("Georgia", size(vlarge) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(2)) ytitle("", size(zero) margin(zero) axis(2)) yscale(lstyle(none) titlegap(0) outergap(0) axis(3)) ytitle("", size(zero) margin(zero) axis(3))) ///
, ///
legend(region(lstyle(none) lcolor(white)) order(4 "Supporters" 3 "Opponents") size(*1) symxsize(*.3) symysize(*.3) forcesize ring(0) position(11) colgap(tiny) row(2)) plotr(m(medsmall)) ///
)

graph save fig5_geo, replace

graph combine fig5_lat.gph fig5_geo.gph, xsize(20) ysize(10) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) b1("Treatments" , size(*1.2) margin(0 0 0 0) color(black)) l1("Support for war" , size(*1.2) margin(0 0 0 0) color(black)) rows(1) cols(3) 


graph export figure_4.png, width(4000) replace
graph export figure_4_pdf.pdf, replace



********************************************************************************
************************************ Figure 5 **********************************
********************************************************************************

* Deployment of missiles *
use Deplyment_of_missiles, clear

twoway (bar freq date, color(gs10%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=19, color(purple%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=23, color(purple%70) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date ==11, color(black%99) mlwidth(thin) lwidth(thin) msize(*.75)) ///
, ///
ylabel(#5, labcolor(black) labels labsize(*.7) angle(horizontal) valuelabel nogrid noticks nogextend) ///
ytitle(" ", margin(zero)) ///
xlabel(1(1)24, labcolor(black) labsize(*.7) angle(vertical) noticks valuelabel) ///
xtitle("") ///
plotr(m(zero)) ///
title({bf:Deployment of missiles}, size(*.7) color(black) margin(1 1 3 1) position(12)) ///
legend(off) ///
graphregion(fcolor(white) lcolor(white) icolor(white) margin(2 2 2 2)) ///
plotregion(fcolor(white) lcolor(white) icolor(white)) ///
xsize(20) ysize(16)


graph save Deplyment_of_missiles, replace


* Threat to Russia *
use Threat_to_Russia, clear

twoway (bar freq date, color(gs10%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=19, color(purple%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=23, color(purple%70) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date ==11, color(black%99) mlwidth(thin) lwidth(thin) msize(*.75)) ///
, ///
ylabel(#5, labcolor(black) labels labsize(*.7) angle(horizontal) valuelabel nogrid noticks nogextend) ///
ytitle(" ", margin(zero)) ///
xlabel(1(1)24, labcolor(black) labsize(*.7) angle(vertical) noticks valuelabel) ///
xtitle("") ///
plotr(m(zero)) ///
title({bf:Threat to Russia}, size(*.7) color(black) margin(1 1 3 1) position(12)) ///
legend(off) ///
graphregion(fcolor(white) lcolor(white) icolor(white) margin(2 2 2 2)) ///
plotregion(fcolor(white) lcolor(white) icolor(white)) ///
xsize(20) ysize(16)

graph save Threat_to_Russia, replace


* Russian security *
use Russian_Security, clear

twoway (bar freq date, color(gs10%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=19, color(purple%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=23, color(purple%70) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date ==11, color(black%99) mlwidth(thin) lwidth(thin) msize(*.75)) ///
, ///
ylabel(#5, labcolor(black) labels labsize(*.7) angle(horizontal) valuelabel nogrid noticks nogextend) ///
ytitle(" ", margin(zero)) ///
xlabel(1(1)24, labcolor(black) labsize(*.7) angle(vertical) noticks valuelabel) ///
xtitle("") ///
plotr(m(zero)) ///
title({bf:Russian security}, size(*.7) color(black) margin(1 1 3 1) position(12)) ///
legend(off) ///
graphregion(fcolor(white) lcolor(white) icolor(white) margin(2 2 2 2)) ///
plotregion(fcolor(white) lcolor(white) icolor(white)) ///
xsize(20) ysize(16)

graph save Russian_Security, replace



* Russian speakers' security *
use Russianspeakers_security, clear

twoway (bar freq date, color(gs10%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=19, color(purple%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=23, color(purple%70) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date ==11, color(black%99) mlwidth(thin) lwidth(thin) msize(*.75)) ///
, ///
ylabel(#5, labcolor(black) labels labsize(*.7) angle(horizontal) valuelabel nogrid noticks nogextend) ///
ytitle(" ", margin(zero)) ///
xlabel(1(1)24, labcolor(black) labsize(*.7) angle(vertical) noticks valuelabel) ///
xtitle("") ///
plotr(m(zero)) ///
title({bf:Russian speakers' security}, size(*.7) color(black) margin(1 1 3 1) position(12)) ///
legend(off) ///
graphregion(fcolor(white) lcolor(white) icolor(white) margin(2 2 2 2)) ///
plotregion(fcolor(white) lcolor(white) icolor(white)) ///
xsize(20) ysize(16)

graph save Russianspeakers_security, replace


* Concern *
use Worry_Concern, clear

twoway (bar freq date, color(gs10%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=19, color(purple%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=23, color(purple%70) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date ==11, color(black%99) mlwidth(thin) lwidth(thin) msize(*.75)) ///
, ///
ylabel(#5, labcolor(black) labels labsize(*.7) angle(horizontal) valuelabel nogrid noticks nogextend) ///
ytitle(" ", margin(zero)) ///
xlabel(1(1)24, labcolor(black) labsize(*.7) angle(vertical) noticks valuelabel) ///
xtitle("") ///
plotr(m(zero)) ///
title({bf:Concern}, size(*.7) color(black) margin(1 1 3 1) position(12)) ///
legend(off) ///
graphregion(fcolor(white) lcolor(white) icolor(white) margin(2 2 2 2)) ///
plotregion(fcolor(white) lcolor(white) icolor(white)) ///
xsize(20) ysize(16)

graph save Worry_Concern, replace



* Russophobia *
use Russophobia, clear
drop in 25/39

twoway (bar freq date, color(gs10%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=19, color(purple%50) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date >=23, color(purple%70) mlwidth(thin) lwidth(thin) msize(*.75)) ///
(bar freq date if date ==11, color(black%99) mlwidth(thin) lwidth(thin) msize(*.75)) ///
, ///
ylabel(#5, labcolor(black) labels labsize(*.7) angle(horizontal) valuelabel nogrid noticks nogextend) ///
ytitle(" ", margin(zero)) ///
xlabel(1(1)24, labcolor(black) labsize(*.7) angle(vertical) noticks valuelabel) ///
xtitle("") ///
plotr(m(zero)) ///
title({bf:Russophobia}, size(*.7) color(black) margin(1 1 3 1) position(12)) ///
legend(off) ///
graphregion(fcolor(white) lcolor(white) icolor(white) margin(2 2 2 2)) ///
plotregion(fcolor(white) lcolor(white) icolor(white)) ///
xsize(20) ysize(16)

graph save russophobia, replace

graph save Threat_to_Russianspeakers, replace

graph combine Deplyment_of_missiles.gph Threat_to_Russia.gph Russian_Security.gph Russianspeakers_security.gph  Worry_Concern.gph russophobia.gph, xsize(20) ysize(10) graphregion(fcolor(white) lcolor(white) ifcolor(white) ilcolor(white)) rows(2) cols(3) 

graph export figure_5.png, width(4000) replace
graph export figure_5_pdf.pdf, replace








